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ABSTRACT 

In the standard Big Bang cosmology the canonical value for the ratio 
of relic neutrinos to CMB photons is 9/11. Within the framework of the 
Standard Model of particle physics there are small corrections, in sum about 
1%, due to slight heating of neutrinos by electron/positron annihilations and 
finite-temperature QED effects. We show that this leads to changes in the 
predicted cosmic microwave background (CMB) anisotropies that might be 
detected by future satellite experiments. NASA's MAP and ESA's PLANCK 
should be able to test the canonical prediction to a precision of 1% or better 
and could confirm these corrections. 



Introduction. Neutrinos are almost as abundant as photons in the Universe 
and contribute almost as much energy density . Under the assumption that 
neutrinos decoupled completely before electrons and positrons annihilated (at 
a time of around lsec), the ratio of the number density of neutrinos to that 
of photons is 



n 7 V 11 

where N v = 3 is the number of neutrino species. Further, because of the 
heating of the photons by e + /e~ annihilations, the ratio of the neutrino 
temperature to the photon temperature is (4/11) 1 / 3 = 0.714. It follows that 
the ratio of the energy density of neutrinos to that of photons is 

PJL = 7 (±) 4/i Ny = 0.681. (2) 



p 7 8 Vll 

It has been pointed out that the assumption that neutrinos decoupled 
completely before e + /e~ annihilations is not precisely valid There is 
now a consensus that the neutrinos share in the heating somewhat, so their 
number and energy density is slightly larger than the canonical values, Eqs. 
dD, §). The increase is equivalent to having slightly more than three neutrino 
species and the canonical ratios. (This is just a heuristic device, of course; the 
actual number of generations is fixed at three.) The change in the effective 
number of neutrino generations is f|, |], H; 0, Hi 

5Nl D = 0.03. (3) 

The first calculations 0, |3], |J of this effect were "one-zone" estimates that 
evolved integrated quantities through the process of neutrino decoupling. 
More refined "multi-zone" calculations tracked many energy bins, assumed 
Boltzmann statistics and made other approximations || ||. The latest re- 
finements have included these small effects as well [0, BJ. (A very recent 
calculation makes no approximations whatever and tracks the neutrino mo- 
mentum distribution over 8 orders of magnitude in momentum |J, and arrives 
at slightly higher value, SNl D = 0.045. Until the discrepancy is understood 
we will stick with the earlier estimates; in any case it is simple to rescale our 
results.) 

There is another effect operating at roughly the same time which acts 
in the same direction; it involves finite-temperature QED corrections to the 
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energy density of the electron, positron and photon portion of the plasma 



due to interactions ||T0| , [TT|. This effect decreases the energy density of the 
e ± 7 plasma. Consequently this reduces the amount of energy converted to 
photons when electrons and positrons annihilate thereby slightly raising the 
ratio of the neutrino to photon energy densities. This QED effect can also 
be expressed as an increase in the number of neutrino species [fT~0| , [TT|| 



6N QED = o.Ol. (4) 

Together, incomplete annihilation and QED finite-temperature correc- 
tions lead to an increase in the neutrino energy density over the canonical 
value by slightly more than 1%, corresponding to 

6N V = 0.04. (5) 

These two corrections were initially considered in the context of Big Bang 
Nucleosynthesis. Their net effect is to increase the predicted 4 He abundance 
by a tiny amount, AYp = +5 x 10~ 4 , which given the present observational 
uncertainties, oy p > 10~ 2 , is undetectable and likely to remain so for quite 
some time. (The quantity Yp denotes the primordial mass fraction of He.) 

On the other hand, the small increase in the neutrino energy density can 
have a significant - and potentially detectable effect - on another remnant 
of the Big Bang - the cosmic microwave background (CMB). In particu- 
lar, the anisotropies in the CMB are very sensitive to the epoch of matter- 
radiation equality, which depends on the neutrino energy density. As we shall 
show, the sensitivity is so great that by itself the additional energy density 
in neutrinos should be detectable by the very precise measurements of the 
CMB anisotropy that will by made the two forthcoming satellite experiments, 
NASA's MAP and ESA's PLANCK Surveyor. However, the situation is com- 
plicated somewhat by the fact that predictions for the CMB anisotropies also 
depend on other cosmological parameters. 

The aim of this paper is to address quantitatively the detectability of 
the small increase in neutrino energy density due to incomplete decoupling 
and finite-temperature QED effects. For definiteness, we assume that the 
primordial density perturbations that seed structure formation and lead to 
CMB anisotropy were set during inflation and allow five other parameters to 
vary. They are: baryon density (flp)', Hubble constant (H ); amplitude of 
primordial perturbations; slope of primordial perturbations (n); and epoch 
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of reionization. We find that (i) if these other parameters are held fixed (e.g., 
because they are determined by other measurements), then detectability is 
a sure thing; and (ii) if the other parameters are allowed to vary, then the 
situation is less promising; but, assuming non-linear effects are not a serious 
contaminant and polarization of the CMB anisotropy is also measured with 
precision |T^], these small corrections should be detectable. 



Probing Neutrino Physics with the CMB. Anisotropies in the CMB are best 
characterized by expanding the temperature field on the sky in terms of 
spherical harmonics: 

oo I 

T(M)=E E a lm Y lm {8,<f>). (6) 

1=0 m=-l 

A given theory, specified by the primordial spectrum of perturbations and 
cosmological parameters, makes predictions about the multipole amplitudes, 
the aim's. The predictions take the form of statements about the distribution 
of the aim's. Inflationary theories typically predict that each of these coef- 
ficients is drawn from a Gaussian distribution; as such, the distribution can 
be defined by its variance. Thus, the fundamental predictions of inflationary 
models are 

Q = {a lm a* lm ). (7) 

Much effort has gone into computing the C/'s over the last few years; 
they can be calculated very accurately once the cosmological parameters are 



chosen ||13|| . Viewed simplistically, the results of a CMB experiment are 
estimates of the C/'s, with errors given by AC/. Then, by minimizing a x 2 
statistic 

the underlying set of unknown cosmological parameters {Aj} can be esti- 
mated. 

Of course, we cannot know in advance the values of C;'s that a given 
experiment will measure; however, by knowing what we expect for the AC/'s, 
we can estimate how large the uncertainties in the parameters should be 
("error forecasting"). To do this, we assume that the measured C/'s will be 
close to the true C/'s. Then, by expanding \ 2 around its minimum at {A* rue } 
we can estimate the precision to which a parameter can be determined (for 
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further discussion of "error forecasting" in parameter estimation, see e.g., 
Refs. H): 



x 2 ({a,}) - x 2 ({at}) + ~ -S- (a, - Ar) (a, - \r) 

Z, UAiUAj ^_^true 

= X 2 ({Af ue }) + C VJ (% - Af uc ) (A, - Af ue ) . (9) 



The second- derivative (Fisher) matrix Cy carries information about how 
quickly \ 2 increases as the parameters move away from their true values. 
Therefore, under some reasonable assumptions |]15||, the uncertainties in the 



parameters are determined by this matrix. We are interested only in the 
parameter^ N v . If all the other cosmological parameters are held fixed, it is 
a simple exercise to show that its variance is given by 

< = (io) 

If all other parameters are allowed to vary, then 



a 
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{C- 1 )n„,n„ (11) 



To proceed we need to specify 



Cosmological Model. For definiteness, we take this to be a Cold 
Dark Matter model with Hubble constant H = 50 km sec -1 Mpc -1 , 
baryon density Qb — 0.08, COBE-normalized spectrum of scale invari- 
ant density perturbations (i.e., power- law index n — 1), no reionization, 
and energy density in cold dark matter particles f2cDM = 1 — = 0.92. 
(We assume the simplest inflationary prediction of Q = 1.) 

Experimental Errors. Instead of tying ourselves to a particular ex- 
periment, we assume that the experimental uncertainty is given by 

AC i = I \fSr Cl 1 - /max . (12) 

I OO / ^> /max 



lr The only other parameter estimation paper we are aware of which considers N v as 
a free parameter is Jungman et al. (ll). Their analysis, performed several years ago, 
only considered I < 1000 (then considered optimistic). Further they did not consider 
polarization. Where it is possible to compare with them, our results agree. 
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The error [2/(2/+ 1)] 1/2 Q is the smallest possible given that each multi- 
pole amplitude a\ m can be sampled only 2/+1 times; it is the irreducible 
sampling or cosmic variance. Equation ( |T2"D is obviously a simplifica- 
tion, but we have found it to be a reasonable approximation to the more 
realistic formula fl6j which also accounts for detector noise. Further, it 



allows us to display our results as a function of Z max , which will give a 
clear sense of what angular scales need to be probed. We use a similar 
formula for polarization (with different / max ). For orientation, MAP is 
characterized by / max ~ 1000 and PLANCK by / max ~ 2500. 

• Model Parameters. We allow for variation in five parameters besides 
the neutrino energy density: overall amplitude of the spectrum of den- 
sity perturbations, epoch of reionization parametrized by the optical 
depth back to last scattering r, H , and n. 

Figure 1 shows the angular power spectrum (Cj vs. I) and how it changes 
as each of the six parameters is varied. Figure 2 shows the same for the polar- 
ization power spectrum. Using these derivatives, we can evaluate the Fisher 
matrix and compute the expected error in N v . Our results are summarized 
in Fig. 3. 

The lower set of curves in figure 3 show the one-sigma errors on N v if 
all the other parameters are known. Even without polarizarion information, 
both PLANCK and MAP will detect the predicted 5N„ ~ 0.04. If the other 
parameters are not well determined by other considerations, even a very high 
resolution temperature anisotropy experiment will have difficulty detecting 
the small predicted increase in N v . However, as Fig. 3 illustrates, with po- 
larization information, the prospects are considerably brighter. 

We should note here that other effects may also mimic a change in the 
effective number of neutrino species- particularly any field or particle with a 
relativistic equation of state at the matter-radiation equality epoch. For 
example, the presence of a random magnetic field will contribute to the 
relativistic energy density of the universe. For a given average field strength 
B eq at the radiation-matter equality epoch, one could misconstrue this as an 
effective change in the number of neutrino species: 




flV**«0.03 -T^H W4r 4 , (13) 



5 



where h§o is the Hubble constant normalized to 50 km sec -1 Mpc -1 . If for 
example we constrain SN^f ag - < .01 this translates to B e<l < 6 gauss (fio^so) 2 - 
Assuming that the magnetic field B oc a -2 , where a is the scale factor, then 
this limit is an order of magnitude greater than the Big Bang Nucleosynthesis 



limit [17]. More importantly, a magnetic field this of this order may be 



measured or ruled out by Faraday rotation of the CMB and perhaps 
other CMB measurements [19, 20fl , thus the magnetic field effect may be 
disentangled from 8N U . 

Concluding Remarks. As our analysis shows, future, high-precision CMB 
anisotropy measurements have the potential to measure the cosmic energy 
density in neutrinos to a precision of 1% or better. Such a measurement 
would have significant implications: 

• If = 3, further evidence for the existence of the tau neutrino. Note, 
the tau neutrino has yet to be directly detected in the laboratory. 

• Determination that U N V = 3" by CMB anisotropy would confirm the 
canonical assumption for the energy density in relativistic particles at 
the epoch of big-bang nucleosynthesis, which is an important input 
parameter for these calculations. 

• Confirmation of the standard cosmology prediction that T v /T^ = (4/11) 1 / 3 
to better than 1%. This would test the physics of e + /e~ annihilation 
and neutrino decoupling in the early Universe. 

• Confirmation of two small physics effects that together increase the 
cosmic neutrino energy density by about 1%. In particular, this would 
be the first evidence for finite-temperature QED corrections and a con- 
straint to the strength of neutrino interactions in the early Universe. 

• If a deviation from the expected N v = 3.04 is found, evidence for 
additional relativistic particle species (or magnetic field) present in the 
early Universe or new physics in the neutrino sector (e.g., neutrino mass 
or decay) ||21|| . This would have significant implications for big-bang 
nucleosynthesis, structure formation in the Universe, and elementary- 
particle physics. 

Realizing the full potential of the CMB as a probe of the cosmic neu- 
trino backgrounds will require precision polarization and anisotropy maps 
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out to multipole number 3000. This seems very ambitious and perhaps even 
unattainable. Nonetheless, the potential payoff discussed here makes the goal 
worth striving for. If we have learned nothing else in the years since COBE, 
we have certainly learned that the experimenters have consistently manage 
to surprise theorists by achieving more than was once thought reasonable. 

Finally, we should acknowledge that there is room to improve upon our 
analysis. For example, we have assumed only six cosmological parameters 
and ignored prior information about them. We have not considered the 
adverse role that "secondary anisotropics" that are generated at late times 
might play. There is clearly room for more work on this important subject. 

The CMB spectra used in this work were generated by CMBFAST O]. 
This work was supported by the DOE and the NASA grant NAG 5-2788 at 
Fermilab and by the DoE grant DE-FG02-91ER40690 at Ohio State. 
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Figure 1: Temperature anisotropies and their derivatives. The top pane 
shows the quadropole-normalized CMB anisotropy spectrum as a function 
of multipole moment /, for our baseline CDM model: N v = 3, Qb — 0.08, 
^cdm = 0.92, H Q = 50km sec -1 Mpc -1 , scale-invariant primordial pertur- 
bations, and no reionization. The lower panes show the derivatives of the 
C/'s with respect to the model parameters 9 i: d\nCi/d\n9i. For the optical 
depth to reionization parameter, dln^i/dr is plotted. 
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Figure 2: Same as Fig. 1, but for the P^'s, the electric field polarization 
anisotropies. 



10 



I I I I I I I I I I I 

-100 200 500 1000 2000 



maximum 1 

Figure 3: One-a sensitivity to SN U , for an experiment cosmic-variance limited 
up to some maximum multipole moment. The horizontal line, 5N U = 0.04, is 
the change in effective number of neutrino families due to neutrino heating 
and the QED effect. The bottom two curves are for the case where all 
cosmological parameters except N v are fixed, while the top curves represent 
the case where all parameters are determined simultaneously. For each group, 
the dashed line shows the results using only temperature anisotropy data, 
while the solid line shows the improvement obtained by including polarization 
data in the analysis. 
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